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Abstract. 

In this paper we report calculations of potential energy curves in the 1.2a.M. < R < 
100a.tt. range at Multireference Configuration Interaction (MRCl) level for doubly 
excited states of the H 2 molecule. We have focused on the Q 2 states which lie between 
the second and third ionization thresholds of molecular ion, i.e., state in 
which lie the H(21) -|- H(21’) dissociation channels. The MRCl approach allowed us 
to successfully identify for the first time the molecular state which dissociates into 
hydrogen atoms at 2s state. Further, Generalized Oscillator Strength as a function of 
transferred momentum for three doubly excited states is also presented. 
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1. Introduction 

One of the first evidences of the antoionized donbly excited states of H 2 was obtained in 
the experiment of Crowe and McConkey in 1973 [1] in which electrons collided with H 2 
at impact energies of 29, 33, 50, 100 and 399 eV. The data revealed the existence of new 
nnknown states besides the gronnd state and the hrst excited state of H^. The 
hrst theoretical calculations of the states were made by Bottcher and Docken in 1974 
[2], who calculated the energy position and widths for dissociative Qi ^fl^ and Q 2 
states for la.u. < R < 10a. m., using the Feshbach projection-operators [3]. More recently 
Sanchez and Martin [1] have done systematic calculations in which they estimated tens 
of doubly excited states symmetries and their autoionization widths for „ 

and in the range of 0 < i? < 5a.u.. They employed the Feshbach projection- 

operator method and a representation of the non resonant continuum, with B-spline 
type functions. In what concerns the dissociation into neutral atoms, perhaps the most 
relevant theoretical work in the recent years was done by Dalgarno and collaborators 
[5] who calculated the doubly excited H 2 states converging to H(n=2) -|- H(n’=2) for 
all the internuclear distances. Using different techniques to solve for each internuclear 
distance region and taking special care of intermediate distances in order to connect the 
different regimes in which calculations were done, they were able to calculate potential 
energy curves from 3 a.u. to 200 a.u. Recently our experimental group [6] observed the 
pair H(2s) -|- H(2s) by measuring neutral fragments in coincidence after H 2 dissociation 
induced by electron impact, conhrming for the hrst time the existence of such states. 
The purpose of the present work is to shed light on the description of the doubly 
excited states that may dissociate onto a pair H(2s) -|- H(2s), since this dissociation 
channel, as discussed in [7j, may be used to probe the spin coherence between the two 
hydrogen atoms in future applications. For this, we have used an approach that mixes 
the standard Cl [8] method with MCSCF [9], known as MRCI [10] in the treatment of 
two-electron diatomic molecules allowing us to rise accurate potential curves at every 
internuclear distance. In addition, since we intend to work with our experimental 
group [6], is desirable to provide it with cross section calculations within the First 
Born Approximation, that would be helpful in the experimental apparatus adjustment. 
For this purpose, we begin by presenting Generalized Oscillator Strength calculation for 
some doubly excited states. Finally, it is important to note that the main feature of 
our calculations is to provide for the hrst time an unambiguous characterization of the 
electronic state that goes, in the dissociation limit, to H(2s) -|- H(2s) fragments. 

We start describing the theoretical methods used for the electronic state energies 
and generalized oscillator strength. Results for three molecular regimes, Franck-Condon 
(FC) (1.2a.u. - 3.0a.u.), intermediate (3a.u. - 12a.u.) and asymptotic range (> 80a. m.). 
We also calculate generalized oscillator strength for those states. At last, conclusion 
and perspectives are pointed out. 
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Atomic units {e=me=h=l) are used throughout, except where indicated. The origin of 
the electronic coordinates is located at the center of mass of the molecule and has its 
z-axis lying along the internuclear axis. 

2.1. Electronics States 

In this study we implement the Huzinaga basis-set m, similar to that used by Borges 
et al na, which consists of cartesian gaussian-type (12s, 6p, 3d, lf)/[9s, 6p, 3d, If] 
and provides a total of 110 molecular orbitals: 29ag, 29au, 227rg, 227r„, 46 and 40. The 
electronic states are obtained within the Born-Oppenheimer approximation, making use 
of a hybrid technique, which is a mixture of the Multi-Conhgurational Self-Consistent 
Field (MCSCF) [9] and the Conhguration Interaction (Cl) [8] methods (also called multi¬ 
reference Cl (MRCI) [lO]). We perform a state averaged MCSCF (SA-MCSCF) over a 
set of 20 states, all of the symmetry. These states are used here for two reasons: 
hrst, they allow us to locate the position in energy of the states of interest, especially 
those of the symmetry which dissociate in H(2s) -|- H(2s); second, we are able to 
optimize the orbitals of interest which generate the Q 2 states which have A = |m^| = 2, 
where ruz is the projection of the electronic angular momentum Lz along the internuclear 
axis. For this purpose, we have used an active space composed of 12 molecular orbitals, 
nag (n = 1-3), nau (n = 1-3), mig (n = 1-2), miu (n = 1-2) which comprise the main 
conhgurations of the electronic states from the Franck-Condon to the dissociative region. 
Since we are interested in the Q 2 states, we have excluded the Icr^ and lcr„ orbitals of 
the conhgurations at the Cl step. However, in order to generate both P(|lscrg) -|- |1 s(T„)) 
and Q = (1 — P) [3] precise spaces, such orbitals are still calculated at the MCSCF step. 
Hence, potential energy curves of all symmetries at the {la.u. < R < 100a.u.) range 
have been obtained. 

2.2. Generalized Oscillator Strength 

In this paper we calculate the Generalized Oscillator Strength (GOS) for some doubly 
excited states in order to evaluate the cross section of excitation of the H 2 molecule, 
in the vertical approximation 0. In this work, we perform the electronic part of the 
GOS, which is given by the following expression: 



( 1 ) 


where AE is the transition energy, K is the module of transferred momentum, gn is 
the degeneracy of the excited electronic state and Rq, the distance of equilibrium. The 
integral in dQ is performed over all possible orientations of the molecule with respect 
to K. The amplitude of scattering within the First Born approximation (FBA) for a 
molecule H 2 is given by: 
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eon{K,RQ,VL) = - / driY2f}n{ri,r2]Ro)\^exp{i\^.ri)]iiQ{ri,r2]Ro) ( 2 ) 

i=i 

In order to calculate the generalized oscillator strength, we have implemented the 
quadrature method of Gauss-Hermite to the electronic coordinates integrals, which 
works hne in a set of Cartesian Gaussian basis. The electronic states were obtained 
as previously mentioned, i.e., MRCI approach. One should note the important fact 
that, in a single determinant approach, the ground state differs from doubly excited 
states for more than one excitation which, for a one electron operator, would result in 
a identically zero transition value (Condon-Slater rule [H]). However, our technique 
takes into account highly correlated wavefunctions, in such a way that the states are 
not described by only one Slater determinant, but a combination of them (Cl), resulting 
in a calculation of transition matrices of two excitations in only one step. 

3. Results and Discussion 

Here we present our results for the potential curves of Q 2 states from the Frank-Condon 
region (1.2a.M. < R < S.O.a.u.) to the dissociative region (f^ 100a. m.) as the Generalized 
Oscillator Strength, within the vertical approximation, using the MRCI [TO] . 

3.1. Franck-Condon range 

The FC region is extremely important, especially in collisional processes, since that is 
the region where the electronic transitions take place. From the theoretical point of 
view, a reliable calculation of excitation cross sections requires that the wavefunctions 
are quite accurate in such region. In order to test the quality of our wavefunction, we 
have calculated the dipole moment of the H 2 molecule and show in Table [T] the results 
for the excited states and respectively. Note that the results present an 

excellent agreement with the exact value calculated by Wolniewicz et al [15], with a 
maximum discrepancy of about 0.6%. 

In Fig. in we show the energy values as a function of the internuclear distance for the 
^Hg and ^H^ symmetries obtained from our technique and compared with 
the results of Sanchez and Martin et al Note that there is an overall good agree¬ 
ment and, in the case of the second root of symmetry, our curve lies below that 
of Sanchez and Martin, which from the variational point of view is an indicator of the 
quality of our results. In Fig. |2]a note that the Cl wavefunctions of the Q 2 in FC 
region are mainly composed of the 2p7r^i; uAtt^ conhgurations, clearly converging to the 
H^ state. Figs [I]&, Hjc andUld show some Q 2 states, namely, symmetries 

running parallel to the ionic hydrogen states (as also observed by Sanchez and Martin 
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Figure 1. Potential Energy Curves as a function of the internuclear distance. Dashed 
lines, Sanchez and Martin [i]; triangles, our results: (a) the three hrst Q 2 
symmetries; (b) two first Q 2 ^E+; (c) Q 2 (d) Q 2 


Table 1. Dipole Transition Moment for some internuclear distances calculated by (a) 
Borges [12], (b) Wolniewicz [15] and (c) This work 


State 

R=1.0ao 

R—1.2ao 

R=1.4ao 

R=1.6ao 

bo 

0 

R—2.0a[ 

B'EJ 

0.7269“ 

0.8359“ 

0.9236“ 

1.025“ 

1.125“ 

1.218“ 


0.7650'’ 

0.8708'’ 

0.9821'’ 

1.096'’ 

1.208'’ 

1.313'’ 


0.7654“ 

0.8707“ 

0.9819“ 

1.095“ 

1.207“ 

1.313“ 

B'lEJ 

0.3262“ 

0.3584“ 

0.3886“ 

0.3966'’ 

0.4148“ 

0.4352“ 

0.4355'’ 

0.4478“ 

0.4388'’ 


0.3394“ 

0.3709“ 

0.3990“ 

0.4210“ 

0.4341“ 

0.4366“ 


3.2. Intermediate Range 

The range 3a. m. < R < 12a. m. can still be considered as intermediate regime, in the 
sense that in this region the potential energy wells of the Q 2 states can still be found. In 
this region, there is still electronic correlation, however smaller. This region is important 
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Figure 2. (a) Weight of the main configurations of IQ 2 as a function of 

internuclear distance: in dotted lines (2p7r„;nA7r„); dashed lines correspond to the 
{2s(7g]n\(jg). In (b), (c) and (d), solid lines are the and Hj ion energy, 
dotted and dashed lines are the potential energy curves as a function of the internuclear 
distance of some electronic doubly excited states ignoring the nuclear interaction 


for the verification of crossings between states, which are of fundamental importance in 
calculations of survival probability in a dissociative regime. The states were calculated 
hrstly with a careful choice of the active space, including the 2s(Tg, 2sc^^l, 2p7r„, 2p'Kg, 
SpTiu, SpTig, 3sag and 3sau orbitals. As mentioned before, the doubly excited states are 
characterized by the exclusion of conhgurations with occupation on the orbitals of type 
IcTg (Qi) or lag and lcr„ (Q 2 )- In the Fig. [3l we present the potential energy curves for 
the and states compared with the data of Vanne et al [5] at the 

3.0a.u. < R < 12a.M. range and, as can be seen, our results show a discrepancy smaller 
than 1%. Here we should emphasize that our basis-set is in a number of 110 molecular 
orbitals, whereas Vanne’s calculation took about 150 functions. Values of electronic 
energy of the singlets and triplets for some internuclear distances also are presented in 
Table [2l It must be noted that the advantage of our technique is that we are able to 
work in all internuclear regions without concerning about the adequacy of the basis-set 
to a certain internuclear region, which is not possible in the case of Hartree-Fock based 
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Table 2. Some values of electronic energy for the Q 2 states (a) This work, (b) Vanne 
et al [5] and (c) Sanchez and Martin [4] 


State 

R=3.0ao 

R=4.0ao 

R=5.0ao 

R=6.0ao 

State 

R=3.0ao 

R=4.0ao 

R=5.0ao 

R=6.0ao 

HE+ 

-0.20900“ 

-0.21167^ 

-0.20694“ 

-0.24715“ 

-0.24969'' 

-0.24595“ 

-0.26463“ 

-0.26679'' 

-0.26350“ 

-0.27781“ 

-0.27929'' 

-0.27491“ 


-0.13314“ 

-0.13404'' 

-0.13201“ 

-0.18855“ 

-0.18940'' 

-0.18717“ 

-0.21811“ 

-0.21896'' 

-0.21584“ 

-0.24277“ 

-0.24371'' 

-0.23600“ 


-0.13883“ 

-0.13888^ 

-0.13618“ 

-0.18658“ 

-0.18686'' 

-0.18469“ 

-0.21412“ 

-0.21508'' 

-0.21250“ 

-0.24170“ 

-0.24338'' 

-0.23579“ 

23S+ 

-0.11065“ 

-0.11086'' 

-0.10778“ 

-0.17072“ 

-0.17110'' 

-0.16762“ 

-0.20604“ 

-0.20702'' 

-0.20193“ 

-0.22782“ 

-0.22860'' 

-0.22502“ 


-0.10606“ 

-0.10958^ 

-0.10816“ 

-0.15542“ 

-0.15715'' 

-0.15465“ 

-0.20582“ 

-0.20679'' 

-0.20082“ 

-0.23386“ 

-0.23536'' 

-0.23150“ 


-0.08383“ 

-0.09039'' 

-0.08971“ 

-0.13517“ 

-0.14035'' 

-0.13960“ 

-0.19006“ 

-0.18971'' 

-0.18343“ 

-0.21393“ 

-0.21329'' 

-0.20770“ 


-0.17898“ 

-0.18051^ 

-0.18817“ 

-0.21929“ 

-0.22099'' 

-0.21884“ 

-0.23739“ 

-0.23896'' 

-0.23629“ 

-0.25124“ 

-0.25327'' 

-0.24904“ 


-0.20372“ 

-0.20540'' 

-0.20225“ 

-0.23918“ 

-0.24084'' 

-0.23739“ 

-0.25745“ 

-0.25885'' 

-0.25323“ 

-0.27550“ 

-0.27663'' 

-0.26958“ 

2in„ 

-0.12364“ 

-0.12739^ 

-0.12497“ 

-0.18319“ 

-0.18740'' 

-0.18456“ 

-0.21884“ 

-0.22278'' 

-0.21922“ 

-0.23519“ 

-0.23782'' 

-0.23442“ 

23n„ 

-0.13281“ 

-0.13397'' 

-0.13171“ 

-0.20133“ 

-0.20239'' 

-0.19904“ 

-0.23889“ 

-0.23994'' 

-0.23766“ 

-0.24837“ 

-0.24944'' 

-0.24772“ 

HA, 

-0.21188“ 

-0.21453^ 

-0.21060“ 

-0.24838“ 

-0.25106'' 

-0.24768“ 

-0.26025“ 

-0.26306'' 

-0.26044“ 

-0.26234“ 

-0.26553'' 

-0.26332“ 


-0.12873“ 

-0.12966'' 

-0.12851“ 

-0.18491“ 

-0.18595'' 

-0.18470“ 

-0.21369“ 

-0.21468'' 

-0.21349“ 

-0.22964“ 

-0.23054'' 

-0.22936“ 


Cl approach [1] . Another advantage over Cl with B-spline functions methods, as we will 
see on the next section, is that by using a relaxed basis we are able to uniquely identify 
the way each state dissociates, even in the degenerate cases. 

3.3. Asymptotic Range 

In this internuclear distance, the electronic correlation is almost zero and only a few 
dispersion terms of the order of 1/i?® may play a role between the atoms via van der 
Waals interaction ng. As said before, one of the fundamental differences between the 
MCSCF and HF-based Cl methods is that in the latter, the orbitals remain hxed while 
only expansion coefhcients of the conhgurations are taking as variational parameters. 
However, in MCSCF the orbitals and Cl coefhcients are both optimized. This relaxation 
of the orbitals allows us, through the analysis of the conhguration coefhcients of higher 
weight, to determine whether the molecule is going or not to dissociate on the H(2s) 
+ H(2s) fragments. The main conhguration of the whose dissociation energy 

is -0.25 a.u. (two hydrogen atoms with n = 2) in the asymptotic region, is composed 
by 2sag and 2scr„ molecular orbitals. The shape of these molecular orbitals is shown in 
Fig. m at R=80 a.u., in which the electronic amplitudes are plotted. 3D plots for the 
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electronic amplitude are shown for both the 2s(jg (on the left) and 2s(Ju (on the right) 
orbitals. It is evident from these pictures that, at the separated atoms regime, both 
atomic orbitals present one node which is characteristic of a 2s orbital. 






Figure 3. Potential Energy Curves as a function of the internuclear distance. Dashed 
lines for Vanne et al curves and triangles for this work: (a) Q 2 symmetries,(b) 

Q 2 symmetries, (c) Q 2 and symmetries and (d) Q 2 ^Ilg and ^n„ 
symmetries 

It is important to note that neither the work of Sanchez and Martin HZ] nor any 
other theoretical work on doubly excited states of H 2 molecule were able to provide 
information about the Q 2 electronic states in the dissociation limit, given that the 
states with symmetry may dissociate in H(2/)+H(2Z'), where 1 and V are either s 
or p, that is, they are degenerated and that requires, from the theoretical perspective, 
more information than simply the energy to be able to determine which fragments were 
produced. 

3.4- Generalized Oscillator Strength 

In Fig. [5] is presented the generalized oscillator strength as function of the transferred 
momentum for three doubly excited states of lowest energies (IQ 2 2 Q 2 and 
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Figure 4. In the left, 3D electronic amplitude plot of 2sUg] in the right, the same for 
2sau orbital. Colour scheme: blue line stands for a positive amplitude value and red 
line for negative amplitude (online version) 


IQ 2 of the H 2 molecule. MRCI wavefunctions have been used and the calculations 
have been done at the distance of equilibrium (1.4 a.u.) of the molecule. Note that, 
there are no previous results to be compared. In the Fig. |5]a, for the IQ 2 it can 
be seen that the GOS curve presents a maximum peak at = O.da.n.. This is a very 
useful information that can be used to set an optimum experimental apparatus since 
the GOS is directly related to the excitation cross-section. For the same reason of the 
early state, Fig. |5]6 shows that the transition for the second is forbidden, presenting 
a smoother peak at = l.Sa.n.. Finally, for the state the transition is allowed. 
The corresponding prohle is shown in Fig. 5c. 

4. Conclusions 

One of our goals, namely, to map the Q 2 state which dissociates in H(2s) + H(2s) 
was achieved successfully. Here, again, we should emphasize that no other technique 
was capable to provide this crucial information in what concerns the applications in 
twin atoms [6]. Further, we are able to plot energy curves up to a higher internuclear 
distance than that of Sanchez and Martin UTl, given that their B-spline functions present 
convergence problems for distances greater than 6 a.u.. Although our base is composed 
by Gaussian functions centered at the nuclei, it was still possible to reach convergence 
at distances higher than 100 a.u.. We have also performed calculations for the transition 
moments for the states and and achieved precise results and agreement 

comparing to previous calculations and the exact value [15]. We also were able to 
perform matrix elements calculations in the FBA for some doubly excited states and 
here it is worth pointing out that, as said before, a transition calculation for a doubly 
excited state in the FBA (which is an one-electron operator), should be identically 
null according to Gondon-Slater rule, for a single determinant approach but, our highly 
correlated wavefunction allowed us (despite of an one-electron operator calculation) to 
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Figure 5. Generalized Oscillator Strength as a function of the transferred momentum 
distance for doubly excited states within MRCI method: (a) IQ 2 ^51+ symmetry, 
(b) 2 Q 2 symmetry and (c) IQ 2 symmetry 


be able to lead to a set of non-null values for the transition matrix. This gives us a solid 
and promising basis for the future cross section calculations as well as other physical 
quantities of interest from the theoretical and experimental point of view. We intend 
in the future to rehne this approach, in order to extract more information about the 
excitation cross-section. 
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